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ABSTRACT 


When  a  circularly  polarized  wave  is  incident  on  a  trihedral  corner 
reflector,  the  backscattered  wave  is  circularly  polarized  with  the  oppo¬ 
site  sense  of  rotation.  However,  by  coating  one  of  the  conducting  sur¬ 
faces  of  the  reflector  with  dielectric  layers,  it  is  possible  to  obtain  a 
circularly  polarized  backscatter  wave  having  the  same  sense  of  rotation 
as  the  incident  wave.  With  the  aid  of  a  digital  computer,  the  thicknesses 
of  the  dielectric  layers  can  be  cptimized  in  an  efficient  manner  to  obtain 
the  best  performance  over  the  pertinent  range  of  incidence  angles.  This 
report  includes  a  computer  program  for  this  purpose  and  describes  an 
optimized  four-layer  reactive  wall. 
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COMPUTER  OPTIMIZATION  OF  A  CORNER  RFFLECTOR  FOR 
ECHO  ENHANCEMENT  WITH  CIRCULAR  POLARIZATION 

1.  INTRODUCTION 

This  report  considers  the  design  of  a  reactive  wall  consisting 
of  a  dielectric  multilayer  on  a  conducting  plane.  This  reactive  wall 
is  to  be  used  as  one  face  of  a  trihedral  corner  reflector  in  which  the 
two  remaining  faces  are  uncoated  conducting  planes.  When  a  right- 
circularly  polarized  wave  is  incident  on  the  corner  reflector,  it  is 
desired  in  some  applications  that  the  backscattered  wave  also  be  right- 
circularly  polarized. 

When  there  are  only  one  or  two  dielectric  layers,  the  appropriate 
thicknesses  of  the  layers  can  be  determined  by  graphical  constructions 
on  the  Smith  chart,  using  a  method  devised  by  Kennaughpj.  However, 
Kennaugh  and  Chung[2]  have  shown  that  even  the  carefully  designed 
reactive  walls  have  rather  poor  characteristics  when  there  are  only 
one,  two  or  three  layers. 

Although  the  Smith  chart  technique  yields  a  suitable  design  for  any 
given  angle  of  incidence,  it  does  not  lead  to  a  design  which  is  optimum 
over  a  range  of  angles  of  incidence.  This  report  describes  a  method 
for  obtaining  such  an  optimum  design  with  the  aid  of  an  automatic  digital 
computer.  The  technique  is  similar  to  that  developed  by  Kotik  and 
Cope[  3]  for  the  optimization  of  multilayer  radome  walls. 

The  following  sections  describe  the  technique  and  present  a 
computer  program  and  typical  results  for  an  optimized  reactive  wall 
having  four  layers. 


II.  COMPUTER  OPTIMIZATION 

Consider  a  harmonic  plane  wave  to  be  incident  on  a  perfectly 
conducting  plane  coated  with  a  dielectric  multilayer  as  shown  in  Fig. 

1.  The  reflection  coefficients  for  the  TE  and  TM  cases  can  be  calcu¬ 
lated  with  the  aid  of  a  computer  program  included  in  Report  1 96 8 -  1  [ 4j  . 
If  the  dielectric  layers  are  lossless,  the  TE  and  TM  reflection  coeffici¬ 
ents  are  given  by 
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Fig.  1.  A  plane  wave  incident  on  a  reactive  wall 
consisting  of  four  dielectric  layers  on  a 
conducting  plane. 


(1)  Re  =  eJb 

(2)  Rm  =  e>b' 


That  is,  total  reflection  occurs  and  the  reflection  phase  angles  are 
denoted  by  6  and  t)'.  Let  Re  denote  the  ratio  of  the  electric  field 
intensities  of  the  reflected  and  incident  plane  waves  in  the  TE  case 
(perpendicular  polarization)  and  let  Rm  denote  the  ratio  of  the  magnetic 
field  intensities  of  the  reflected  and  incident  plane  waves  in  the  TM  case 
(parallel  polarization).  If  the  incident  wave  has  right-circular  polariza¬ 
tion,  the  reflected  wave  will  in  general  have  both  right-  and  left-cir- 
cularly  polarized  components.  In  the  present  application  we  are  inter¬ 
ested  in  the  "same-sense  power  reflection  coefficient"  defined  by  the 
ratio  of  the  power  density  of  the  right-circularly  polarized  component 
of  the  reflected  wave  and  the  power  density  of  the  right-circularly 
polarized  incident  wave.  This  same- sense  power  reflection  coefficient, 
denoted  by  R,  is  given  by 

(3)  R  =  cos*  <p 
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where 


(4)  9  =  (6  -  b')/Z 


It  may  be  noted  that  the  reflection  phase  difference  (b  -  &')  is  equal  to 
the  phase  angle  of  the  ratio  of  the  complex  reflection  coefficients: 


(5)  b  -  5'  =  Phase  (Re/Rm)  • 

If  the  dielectric  layers  are  dissipative,  the  same-sense  power  reflection 
coefficient  is  given  by 

(6)  R  =  0.25  |Re  +  Rm | 2  . 


Suppose  that  the  dielectric  constants  «j,  c2,  «j,  and  c«  of  the  four 
layers  are  specified  and  we  are  to  design  the  thicknesses  dlt  d2,  dj , 
and  d*  to  maximize  the  same- sense  power  reflection  coefficient  R  over 
a  given  range  of  incidence  angles.  (In  the  case  of  the  trihedral  corner 
reflector,  Kennaugh[l]  has  shown  that  the  pertinent  range  of  incidence 
angles  is  from  6  =  35°  to  0  =  75°.  )  To  be  more  specific,  suppose  that 
the  layer  thicknesses  are  to  be  designed  to  maximize  the  minimum  value 
of  R  over  this  range  of  incidence  angles.  That  is,  the  "figure  of  merit" 
of  any  given  design  is  taken  to  be  the  minimum  value  of  R  in  the  range 
from  Q  =  35°  to  75°,  and  this  quantity  is  denoted  by  Rm.  Rm  is  a  function 
of  the  thicknesses  dj,  d2>  dj,  and  cU.  It  is  also  a  function  of  the  frequency 
and  the  dielectric  constants  of  the  layers,  but  these  are  assumed  to  be 
specified  and  the  permeability  of  each  layer  is  taken  to  be  that  of  free 
space. 


The  most  straightforward  (but  impractical)  method  for  determining 
the  layer  thicknesses  is  to  program  the  computer  to  calculate  and  print 
out  R  versus  incidence  angle  for  a  large  number  of  designs,  covering  a 
given  range  in  each  of  the  thicknesses.  A  study  of  the  resulting  data  would 
indicate  the  region  of  the  best  design.  However,  this  simple  method  re¬ 
quires  the  calculation  and  inspection  of  the  results  for  ten  thousand  dif¬ 
ferent  designs  obtained  by  taking  ten  increments  in  each  of  the  four  thick¬ 
nesses.  A  more  sophisticated  approach  can  lead  to  the  optimum  design 
with  much  greater  efficiency. 
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This  digital  optimization  technique  proceeds  as  follows.  An 
initial  design  is  specified  in  the  form  of  a  set  of  numerical  values  for 
the  four  thicknesses.  The  computer  is  programmed  to  calculate  the 
figure  of  merit  Rm  for  this  initial  design  and  for  four  slightly  different 
designs  involving  small  changes  in  the  layer  thicknesses.  Based  on 
this  information,  the  computer  determines  the  partial  derivatives  of 
Rm  with  respect  to  the  thickness  of  each  layer.  It  then  arrives  at  a 
new  and  better  design  (having  a  larger  value  of  Rm)  by  changing  the 
thickness  or.  each  layer  by  an  amount  proportional  to  the  partial  deri¬ 
vative  of  Rm  with  respect  to  that  thickness.  For  this  new  design,  Rm 
and  the  partial  derivatives  are  again  calculated  in  order  to  arrive  at  a 
still  better  design.  After  repeating  this  procedure  several  times,  the 
computer  soon  arrives  at  the  optimum  design.  In  a  typical  case,  the 
computer  investigates  40  different  designs  in  order  to  find  the 
optimum.  In  finding  the  partial  derivatives,  it  investigates  four  other 
designs  in  the  vicinity  of  each  of  these  40  basic  ones,  making  a  total 
of  200.  If  this  is  compared  with  the  10,  000  designs  which  must  be 
investigated  in  the  method  described  previously,  the  tremendous  ad¬ 
vantage  of  the  computer  optimization  technique  becomes  apparent. 

The  computer  programs  for  the  optimization  process  are  presented 
in  the  Appendix. 


III.  AN  EXAMPLE  OF  REACTIVE  WALL 

OPTIMIZATION 

As  an  example,  consider  the  optimization  of  a  four-layer  reactive 
wall  having  lossless  layers  with  relative  dielectric  constants  given  by 
«i  =  1,  «2  =  9.6,  <3  =  1,  and  C4  =  9.6.  The  initial  design  is  specified  by 
the  thicknesses  di  =  0.  05\  ,  d2  =  0.  1\  ,  d3  =  0.  05X  ,  and  d*  =  0.  IX.  where 
X  denotes  the  free-space  wavelength.  The  digital  computer  calculates 
the  same-sense  power  reflection  coefficient  R  for  this  initial  design  at 
angles  of  incidence  from  0  =  36*  to  7  5°  in  increments  of  6*.  It  also  cal¬ 
culates  R  versus  0  for  four  slightly  different  designs  in  which  the  thick¬ 
ness  of  one  layer  is  increased  by  0.  001X  .  Table  I  lists  the  thicknesses 
of  the  layers  for  each  of  these  five  cases. 


*  The  initial  design  may  be  based  on  experience  or  on  the  results  of 
a  graphical  solution  on  the  Smith  chart.  To  a  certain  extent,  the 
final  results  will  be  the  same  regardless  of  the  starting  point,  but 
the  computational  costs  are  reduced  by  selecting  a  starting  point 
close  to  the  optimum  design. 
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TABLE  I 

LIST  OF  LAYER  THICKNESSES 
d/k 


Layer  No. 

Case  1 

Case  2 

Case  3 

Case  4 

Case  5 

1 

0.  050 

0.  051 

0.  050 

0.  050 

0.  050 

2 

0.  100 

0.  100 

0.  101 

0.  100 

0.  100 

3 

0.  050 

0.  050 

0.  050 

0.  051 

0.  050 

4 

0.  100 

0.  100 

0.  100 

0.  100 

0,  101 

Table  II  lists  the  power  reflection  coefficient  R  versus  angle  of  incidence 
for  these  five  cases. 


TABLE  II 

SAME-SENSE  POWER  REFLECTION  COEFFICIENT 


9  Case  1  Case  2 


35° 

.3600 

.3425 

40 

.7516 

.7297 

45 

.  9948 

.9913 

50 

.9316 

.  9374 

55 

.7947 

.  7987 

60 

.6913 

.  6925 

65 

.6320 

.6313 

70 

.6150 

.6134 

75 

.6467 

.  6448 

Case  3  Case  4  Case  5 


.  3256 

.3456 

.  3378 

.7045 

.7361 

.7218 

.  9850 

.  993  0 

.  9890 

.9479 

.  9332 

.  9428 

.8120 

.7935 

.  8076 

.7051 

.6881 

.7030 

.6436 

.  6282 

.  6433 

.  6261 

.6115 

.6269 

.6583 

.6441 

.6596 

In  each  case  the  minimum  reflection  coefficient  occurs  at  9  =  35*.  The 
figure  of  merit  for  each  case  is  Rm  =  0,  3600,  0.  3425,  0.  3256,  0.  3456, 
and  0.  3378.  The  best  design  out  of  these  five  cases  is  Case  I,  since  it 
has  the  highest  figure  of  merit.  Since  Rm  decreases  when  the  thickness 
of  any  layer  is  increased  slightly,  it  is  obvious  that  a  better  design  can 
be  obtained  by  decreasing  the  thicknesses  of  all  the  layers,  unless  by 
chance  the  initial  design  (Case  1)  is  optimum. 

From  Table  II  it  is  seen  that  a  change  in  the  thickness  of  layer 
2  has  the  greatest  effect  on  Rm»  and  a  change  in  layer  3  has  the  least 
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effect.  Thus,  in  arriving  at  a  better  design,  we  should  decrease  d2  by 
a  relatively  large  amount  and  layer  3  by  a  relatively  small  amount. 

If  the  figures  of  merit  for  the  five  cases  are  denoted  by  Rj,  R2, 
R*  *  R**  and  R5  ,  the  partial  derivatives  of  Rm  with  respect  to  the  layer 
thicknesses  are  given  approximately  by 


(7) 

9Rm 

3dx 

=  (R2  -Ri)/s 

(8) 

9Rm 

=  (Rj  -  Rx)/s 

8d2 

and  so  forth,  where  s  represents  the  increment  in  thickness  which  is 
0.  001X  in  this  example.  In  order  to  obtain  an  improved  design,  all 
four  layers  in  the  initial  design  are  now  changed  in  thickness  by  amounts 
proportional  to  the  partial  derivatives.  In  this  manner  the  computer 
arrives  at  the  first  improved  design  with  thicknesses  d*  =  0.  0498X  , 
d2  =  0.  0996X  ,  dj  =  0.  0498X  ,  and  d4  =  0.  0997X  .  The  minimum  reflection 
coefficient  over  the  specified  range  of  incidence  angles  is  found  to  be 
Rm  =  0.  3667.  This  represents  a  definite  improvement  over  the  initial 
design,  since  the  figure  of  merit  has  increased  from  0.  3600  to  0.3867. 

After  40  cycles  of  this  procedure,  the  computer  arrived  at  a 
design  which  is  very  close  to  the  optimum.  This  optimum  design  is 
specified  in  Fig.  2  which  also  show's  the  same-sense  power  reflection 
coefficient  as  a  function  of  the  angle  of  incidence.  This  reflection  co¬ 
efficient  exceeds  0.  965  over  the  specified  range  of  incidence  angles. 

Thus,  the  figure  of  merit  of  the  final  design  is  0.  965.  The  entire 
optimization  procedure  required  only  3  minutes  on  an  IBM  7094  com¬ 
puter,  including  a  second  run  in  which  the  thickness  increment  was 
reduced  to  0.  0002X  and  the  incidence  angle  increment  was  reduced  to 
2°.  The  initial  design  for  the  second  run  was  taken  to  be  the  best  de¬ 
sign  obtained  from  the  first  run. 

Figure  2  also  shows  the  same- sense  power  reflection  coefficient 
for  the  mori  practical  case  where  the  losses  in  the  dielectric  layers  are 
taken  into  account.  The  loss  tangent  was  taken  to  be  zero  for  layers  1 
and  3  and  0.  001  for  layers  2  and  4.  The  optimum  layer  thicknesses  were 
found  to  be  0.  0426X  ,  0.  0897X  ,  0.  0464X  ,  and  0.  0980X  in  this  <  ase,  and 
the  minimum  power  reflection  coefficient  is  0.  929.  These  thicknesses 
differ  very  little  from  those  of  the  optimized  lossless  multilayer  which 
are  listed  in  Fig.  2. 


6 


SAME-SENSE  POWER  REFLECTION  COEFFICIENT 


t 


35  40  45  50  55  60  65  70  75 

ANGLE  OF  INCIDENCE  9  (DEGREES) 


Fig.  2.  Same-sense  power  reflection  coefficient  versus 
angle  of  incidence  for  a  circularly-polarized 
wave  incident  on  a  four-layer  reactive  wall 
with  optimized  thicknesses. 


The  optimization  procedure  described  above  was  carried  out 
for  three  reactive  wall  designs  using  lossless  layers  with  different 
dielectric  constants.  The  dielectric  constants  and  the  optimum  layer 
thicknesses  for  these  cases  are  listed  in  Table  III. 
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TABLE  III 

OPTIMIZED  REACTIVE  WALL  DESIGNS 
V/ITH  LOSSLESS  LAYERS 
Design  Frequency:  33.2  Gc 


Dielectric  Constants 


Layer  No. 

Case  1 

Case  2 

Case  3 

1 

1.0 

1.0 

1.0 

2 

6.0 

8.0 

9.6 

3 

1.0 

1.0 

1.0 

4 

6.0 

8.0 

9.6 

Thicknesses 

in  Inches 

Layer  No. 

Case  1 

Case  2 

Case  3 

1 

0.  02752 

0.01810 

0.  01512 

2 

0.  04157 

0. 03547 

0.  03182 

3 

0.  02652 

0.  01860 

0.0164b 

4 

0.  03110 

0.  03578 

0. 03477 

Each  of  these  three  designs  was  found  to  have  excellent  characteristics 
in  the  frequency  range  from  33.  0  to  33.4  Gc  and  the  incidence  angle  range 
from  35*  to  75*.  Of  these  three  designs,  Case  1  is  probably  the  most 
promising  because  the  allowable  tolerances  in  the  la^er  thicknesses  are 
undoubtedly  less  critical  than  in  the  other  two  cases.  Table  IV  lists  the 
reflection  coefficients  for  Case  1  as  a  function  of  frequency  and  angle  of 
incidence. 
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TABLE  IV 

SAME-SENSE  POWER  REFLECTION  COEFFICIENT 

FOR  CASE  1 


Angle  of 
Incidence 

33.  0  Gc 

33.  1  Gc 

33.2  Gc 

33.3  Gc 

33.  4  Gc 

35* 

.81 

.81 

.81 

.81 

.80 

37 

.  88 

.  88 

.89 

.89 

.89 

39 

.93 

.93 

.94 

.94 

.94 

41 

.96 

.97 

.97 

.97 

.98 

43 

.98 

.98 

.99 

.99 

.99 

45 

.99 

.99 

.99 

.99 

.99 

47 

.99 

.99 

.99 

.99 

.99 

49 

.99 

.99 

.99 

.99 

.99 

51 

.99 

.99 

.99 

.99 

.99 

53 

.99 

.99 

.99 

.99 

.99 

55 

.99 

.99 

.99 

.99 

.99 

57 

.99 

.99 

.99 

.99 

.99 

59 

.99 

.99 

.99 

.99 

.99 

61 

.98 

.99 

.99 

.99 

.99 

63 

.98 

.99 

.99 

.99 

65 

.94 

.97 

.99 

.99 

.99 

67 

.91 

.95 

.98 

.99 

.99 

69 

.86 

.92 

.96 

.99 

.99 

71 

.81 

.87 

.93 

.97 

.99 

73 

.74 

.82 

.  88 

.94 

.98 

75 

.67 

.75 

.83 

.90 

.96 

IV.  CONCLUSIONS 

In  a  reactive  wall  for  a  corner  reflector,  the  thicknesses  of  the 
diele  ctnc  layers  can  be  designed  most  efficiently  with  a  computer 
optimization  technique.  Starting  with  a  specified  initial  design,  an 
automatic  digital  computer  is  programmed  to  take  small  increments 
in  the  thicknesses  in  such  a  manner  as  to  move  in  the  direction  of 
maximum  rate  of  improvement.  Excellent  new  designs  are  obtained 
in  this  manner  at  low  cost.  An  optimized  four-layer  wall  is  described 
which  has  a  same-sense  power  reflection  coefficient  exceeding  0.965 
for  all  angles  of  incidence  in  the  range  of  interest. 


9 


Future  research  along  these  lines  should  investigate  the  effects 
of  manufacturing  tolerances  in  the  thicknesses  of  the  various  layers. 

An  improved  d.  sign  should  be  attempted  which  will  have  a  larger  band¬ 
width.  Experimental  verification  should  be  carried  out  for  one  or  more 
of  the  optimized  comer  reflectors. 
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APPENDIX 

DIGITAL  COMPUTER  PROGRAM 


Figure  3  shows  a  computer  program  which  has  been  used  success¬ 
fully  for  optimizing  the  thicknesses  of  the  dielectric  layers  in  a  reactive 
wall.  The  program  is  written  in  the  computer  language  known  as  "Scat- 
ran  .  In  its  present  form  the  program  is  set  up  for  a  four-layer  wall 
with  lossless  layers,  but  it  can  readily  be  modified  for  any  number  of 
lossless  or  dissipative  layers. 

The  input  data  symbols  are  defined  as  follows: 

KKK  =  Number  of  angles  of  incidence 

MMM  =  Maximum  number  of  steps  to  take  in  seeking 

the  optimum  design 

E(I)  =  Relative  dielectric  constant  of  layer  i 

DA,  DB,  DC,  DD  =  Initial  thickness,  in  free-space  wavelengths, 

of  layers  1,  2,  3,  and  4 

DDD  =  Thickness  increment,  in  free-space  wavelengths 

THETA  =  Initial  angle  of  incidence,  degrees 

DTHET  =  Increment  in  angle  of  incidence,  degrees 

The  core  of  the  program,  between  statements  S10  and  SI 00,  is  a 
specialized  version  of  a  program  described  in  Report  1968-1  [4]  for  the 
reflection  coefficients  of  a  plane  dielectric  multilayer.  In  order  to  in¬ 
crease  the  speed  and  the  efficiency  of  the  calculations,  the  general  pro¬ 
gram  given  in  Reference  4  is  specialized  here  for  the  case  where  the 
dielectric  layers  are  lossless  and  have  the  same  permeability  as  free 
space.  Furthermore,  the  general  program  gives  solutions  for  the 
reflection  and  transmission  coefficients  of  a  multilayer  with  and  without 
a  conducting  plane.  The  program  is  simplified  considerably  here  where 
the  multilayer  without  a  conducting  plane  is  of  no  interest. 

The  program  includes  provisions  for  terminating  the  run  if  the 
thickness  of  any  layer  goes  negative.  The  run  also  terminates  if  the 
number  of  "unsuccessful  steps"  exceeds  five.  A  step  is  considered  to 
be  unsuccessful  if  simultaneous  increments  in  all  four  thicknesses  leads 
to  a  new  design  which  has  a  smaller  figure  of  merit  than  that  of  the  pre¬ 
vious  design. 

Just  before  statement  Si 00,  the  quantity  DEL  denotes  the  reflection 
phase  difference  (§  -  &')  which  appears  in  Eqs.  (4)  and  (5),  and  RRR(KK) 
is  the  same- sense  power  reflection  coefficient  which  is  denoted  by  R  in 
Eq,  (3).  Just  after  statement  SI 00,  the  program  directs  the  computer  to 
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prin.‘  out  the  list  of  layer  thicknesses  (JL>1,  D2,  D3,  and  D4)  and  the  power 
reflection  coefficient  RRR(I)  versus  angle  of  ;ncidence. 

The  quantities  Rl,  R2,  R3,  R4,  and  R5  which  appear  below  state¬ 
ment  SI  10  represent  the  figures  of  merit  denoted  by  R*,  R2,  Rj,  etc.  in 
Eqs.  (7)  and  (8)  for  Cases  l  through  5.  These  figures  of  merit  are  used 
to  calculate  the  partial  derivatives  with  -espect  to  the  thickness  of  each 
layer.  Just  above  statement  SI 50,  the  thickness  of  each  layer  is  incre¬ 
mented  by  an  amount  proportional  to  the  partial  derivative  of  the  figure 
of  merit.  The  quantity  D  which  is  calculated  between  statements  SI  10 
and  S150  is  defined  by 


(9)  D  =  |  R2  -  Rl  |  +  I R3  -  Rl  |  +  |  R4  -  Rl  |  +  |  R5  -  Rl !  . 


If  it  were  important  to  take  steps  of  uniform  length  in  moving  toward  the 
optimum  design,  it  would  be  correct  to  define  D  by 


(10)  D=J(R2^  Rlf  +  (R3  -  Rl)2  +  (R4  -  Rl)2  +  (R5  -  Rl)2  . 


However,  the  definition  used  in  Eq.  (9)  and  in  the  computer  program  has 
proven  to  be  satisfactory. 

At  the  end  of  the  computer  program  in  Fig.  3,  a  typical  set  of 
input  data  are  shown. 

Figure  4  shows  the  computer  program  used  for  optimizing  a  four- 
layer  reactive  wall  with  dissipative  layers.  The  list  of  loss  tangents, 
denoted  by  TD(1),  TD(2),  TD(3),  and  TD(4),  must  be  given  along  with 
the  other  input  data  described  previously  for  the  program  in  Fig.  3. 
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***  run 

***  SCATRAN 

C0MPLEX(C1 «  AEC ♦ BEC  .  AMC • BMC  4  EXP i ,EXPI  1 ♦  EXP ,  AECP  «  REC ♦ AMCP  « 
RMC.RSR)- 

D 1 MENS I  ON (E ( 1 00 ) , Z ( 1 00 ) , G C 1 00 ) ,  D 1  < 1 0 ) . D2  C 1 0 ) . D3  ( 1 0 ) . D4 ( 1 0 ) , 
RR( 300 ) ,R(30 )  ,RRR( 30 )  )- 

FIRST  READ  INPUT  ,8, (KKK.MMy  )  — 

N  =  4- 
NN=N+ 1  — 

READ  INPUT, 7 « (E(1).E(2).E(3),E(4))- 
READ  I NPUT , 7 , ( DA , DB • DC , DD , ODD ) - 
READ  INPUT, 7,  (THETA, DTHET  )- 
TPI =6. 2831853- 
E(NN)=1 •- 
RR  <  0 )  =  , 0- 

WRI TE  OUTPUT, 2 « (E(l ),E(2)«E(3)»E(4),E<5) « THETA )- 
WRITE  OUTPUT, 2- 
L  =  1  - 

D 1 ( 1 ) =DA- 
D2( 1 ) =D0- 
D3< 1 )=DC- 
D4( 1 ) =DD- 

DO  T HROUGH ( S 1 50 ) ,MM=1 ,  1  ,MM.LE*MMM- 
DO  THROUGH ( S20 ) , I =2 , 1 , 1 ,LE • 5- 
D 1 (  1  ) =D I  ( 1  )- 
D2(  I  )=D2< 1  )- 
D3< I  >=D3( 1  )- 
S20  D4( I >=D4 ( 1 )- 

D I (2 )=D1 ( 1 ) +DDD- 
D2 ( 3  )  SD2 ( 1 ) +DDD- 
D3 ( 4 ) =D3  < 1  J+DDD- 
D4 (5 ) =D4 ( 1 ) +DDD- 

PR0VIDED(D1 (  1  )  .L.O.  ), TRANSFER  TO  (S160)- 
PR0VIDED(D2( I )«L.O. ), TRANSFER  TO  (SI60)- 
PR0VIDED(D3( 1  )  .L.O.  ), TRANSFER  TO  <S160>- 
PR0VIDED<D4 (  1  )  ,L.O.  ) .TRANSFER  TO  JS160)- 
DO  THROUGH ( S 1  1 0  )  . LL  = 1 ,  1 *LL .LE • 5- 
Z ( 1  )  =TP I *D 1  <LL  )- 
Z ( 2 ) =Z ( 1  )  +  TP  I *D2  <  LL ) - 
Z  <3>=Z<2 )  +  TPl*D3(LL )- 
Z(4  )=Z(3)  +  TPI*D4(LL )- 
THE  T  =  THET  A- 

DO  THROUGH! SI  00 ) *KK=I , 1 , KK .LE • KKK- 
TH= .01 745329*  THET  — 

SS=SIN. (TH)*SIN. (TH)— 

CC=COS. (TH)- 

DO  THROUGH (Si  0 ),.  =  !,  1,1 .LE.N- 
SlO  G(I )=SQRT, (E( I)-SS)- 

G ( NN  )  =  CC- 
AEC  = 1 .- 
BEC  =- 1 •- 
AMC= 1 •- 
9YC= 1 •- 

DO  THROUGH < S50 ) , I  1 =2 . 1  ,  I  1  .LE.NN- 

Fig.  3.  Digital  computer  program  for  optimizing  the  thicknesses 
of  a  reactive  wall  with  four  lossless  layers. 


S50 


SI  GO 


SI  02 
SI  10 


S  I  50 

S  1  60 

#«* 

9  50 

I  .0 
•  05 
35.0 


I  -  I  I  —  I  — 

F=G( I )/G< I  I  )- 
Y I *G ( I  ) *Z ( I  1- 
Y I  I =0 (II  ) #Z ( I  1- 
FP=  1  •  -f  F  — 

FM  = 1 .-F- 

EXPI=COS. (YI ) +. I . S I N. ( Y I  1  — 

EXP I  I =C0S. (YI I  )+. ! .SIN. (YI I  )- 
AECP=AEC- 

amcp=amc- 

AECs (,5/EXPI  I  )»(AECP*EXPI*FP  +  BEC*FM/EXPI  )- 
BEC=.5*£XPI  1*  (AECP*EXPI  *FM-»-BEC*FP/EXPI  )- 
F=F*E ( I  I  )/E( I  1- 
FP= 1 • +F- 
FM= 1 • — F - 

AMC=(.5/EXPI  I  )*(AMCP*EXPI*FP  +  BMC»FM/EXPI  )- 

BMC=.5*£XPI  I  *  ( AMCP*EXPI  *FM-fBMC*FP/EXPI  1- 

RSR=3EC*AMC/( AEC*BMC )- 

DEL=FATAN2* ( • IMAG.RSR. .REAL.RSR)- 

RRR ( KK )  =  • 5* ( 1 •♦COS* (DEL) 1- 

thet=thet+dthet- 

WRITE  OUTPUT. 2* (D1  < LL  1 ♦ 02 (LL ) . 03 ( LL  )  ♦ 04 ( LL ) ) - 
WRITE  OUTPUT. 2. ( ( RRR ( I  ) ♦ I =1 . 1 . I .LE.KKK  )  1- 
R(LL)=RRR(1 )- 

DO  THROUGH (SI  02 )«  I =2 « 1 . I.LE.KKK- 

PROV I  DEO ( RRR ( 1  ) «L  «R ( LL  )  )  ♦ R (LL ) =RRR ( 1  1- 

CONT INUE- 

R1=R( I ) - 

R2=R<2 )- 

R3=R  <3  )- 

R4=R (4  )- 

R5=R (5 ) - 

0= • ASS* (R2-R1 1  +  • ABS • (R3-R1 l+.ABS. (R4-R1 )+*ABS* (R5-R1 1- 

RR(MM)=R1- 

WRITE  OUTPUT .2- 

wr  i  te  output  . 2 .  ( oi  ( 1  ) .  D2  ( 1  i  ♦  03  ( 1 ) .  04  ( 1  1 1  - 

WRITE  OUTPUT .2. (RI .R2.R3.R4. R5 ) - 

write  Output . 2- 

PROV I OED ( R 1 *L »RR ( MM- 1  )  )  .L=L  +  1- 

PROV I  DEO ( L • G • 5 ) »  TRANSFER  T0(S1 60  )- 

01(11=01(1 l+DDD* (R2-R1 1/0- 

02(1) =02  < I l+DDD* (R3-R1  1/0- 

03(1 1=03(1 1 +000* ( R4-R 1 1/0- 

D4 ( 1 1 =D4 ( 1 1 +000* (R5-R1 1/0- 

wr i te  output. 2- 

CONT INUE- 

END  PROGRAM (FIRST  1  - 
DATA 

9,6  WO  9.6 

.10  .05  .10  .001 

5.0 


Fig.  3.  Digital  computer  program  for  optimizing  the  thicknesses 
of  a  reactive  wall  with  four  lossless  layers,  (cont). 
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*#*  RUN 

***  SCATRAN 

COMPLEX (C ! .AEC.9EC.AMC.SMt.EXPI «EXPl 1 ,EXP. AECP.REC « AMCP. 
RMC«EC«G.CSORTL..CEXPL# .F.yI •  Y  I  I *  FP»  FM )  — 

D 1 MENS I  ON (E  < l O0),Z( 1 00  >«G( 100) «D1 ( 10 )«02( 10) t03< 10 ) «D4 ( 10) • 
RR(300) »R  ( 30 ) .RRR(30).TD( 10), EC  I 10)  )- 
FIRST  READ  I NPUT • 3 •  HOCK . MMM ) - 

N=4- 
NN=N+1- 

READ  I NPUT. 7. (EC1).E(2).E(3).E(4)>- 
READ  I NPUT. 7. < TD < 1  ) • TD  < 2 ) . TD <  3 ) . TD < 4 ) ) - 
READ  INPUT. 7. < DA . DB • DC . DD . DDD ) - 
READ  I NPUT • 7 • ( THETA • DTHET ) — 

TP1 =6*2831853- 
C1  =  .  I  •  1 •- 
E(NN)=1 . - 
EC (NN ) = 1 • - 
RR ( 0 ) = . 0- 

WRITE  OUTPUT, 2* (E(1).E(2).E(3).E(4).E(5) .THETA  )- 
WRITE  OUTPUT. 2. < TD ( 1 ) • TD < 2 ) . TD ( 3 ) » TD (4 ) ) - 
WRITE  OUTPUT. 2- 
L  =  1  - 

D 1  ( 1  )  =DA- 
D2( 1 ) =DB- 
D3( 1 ) =DC- 
D4( 1 ) =DD- 

DO  THROUGH (S5). 1=1 . 1 . 1 .LE.N- 
S5  EC( 1  )=E  (  I  )*( 1 .-C1*TD( 1  ) )- 

DO  THROUGH (SI  50 ) • MM= i • 1 .MM.LE.MMM- 
DO  THROUGH ( S20 ) . 1 =2 *  1 . I *LE»5- 
D1  < 1  )=Dl  <  1  )- 
D2 ( I ) SD2 ( 1 >- 
D3( I ) -D3 ( 1 )- 
S20  D4 (  I  )  =D4 ( 1  )  — 

D 1 ( 2 ) =D 1 ( 1 )+DDD- 
D2 ( 3 ) =D2 { 1 ) +DDD- 
D3  ( 4  )  SD3  (  1  )+DDD- 
D4 ( 5 ) =D4 ( 1 )+DDD- 

PROVIDEDIDI ( 1  ).L#0. ) .TRANSFER  TO  <Sl60>- 
PR0VIDED<D2« 1  ).L*0. ) .TRANSFER  TO  (S160)- 
PR0vIDED(D3( 1  ).L.O.  ), TRANSFER  TO  <S160>- 
PR0VIDED<D4 ( 1 )*L.O. ) .TRANSFER  TO  <S160>- 
DO  THROUGH (Si  10) »LL  =  1 . 1 .LL.LE.5- 
Z ( 1  )=TPI#D1 (LL  )- 
Z ( 2 ) =Z ( 1  )  +  TPI*D2<LL  )- 
Z  (3  )  =  Z  <2  )+TPl*D3<l.L  >- 
Z<4 )-Z<3)+TPI*D4(LL)- 
THET=THETA- 

DO  THROUGH! SI  00 ) »KK  = 1 ♦  1  .KK.LE.KKK- 
TH= .0 1 74 5329* THE T- 
SS3  SIN. <TH)#SIN. <TH)- 
CC=COS. (TH)- 

DO  THROUGH(SlO) . 1=1 ♦ 1 . I .LE.N- 
SlO  G( 1 >=CI*CSQRTL. (EC( I )-SS)- 

G(NN)3CI*CC- 
AEC  = 1 .- 
BEC  =  - 1 •- 

Fig.  4.  Digital  computer  program  for  optimizing  the  thicknesses 
of  a  reactive  wall  with  four  lossy  layers. 
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S50 


SI  00 


SI  02 
SI  10 


SI  50 
S  1 60 

**• 

9  20 

1.0 
.000 
•042525 
35.0 


AMC* 1 •— 

BMC*  I  •  — 

DO  THROUGH  <  S50 )  «  1  1  =fe  «  1  «  I  I .lE.NN- 
1=11-1- 
F=G(  I  )/G  (  m- 
VI =G( I )*Z( 1 >- 
VI I«G( I  1 )*Z(  I  )- 
FP= 1 .+F- 
FM* 1 . — F- 

EXP 1 =CEXPL. { Y 1 ) — 

EXP I I»CEXPL. ( Yl  1  >- 
AECP* AEC- 
AMCP= AMC— 

AEC * ( • 5/EXP  1  I  )* <AECP*EXPI*FP+BEC*FM/EXP!  )- 
BEC* .5*EXP I  I* ( A£CP*EXPI*FM+BEC*FP/EXP1  )- 
FsF*E { 1  I  )/E ( 1  >- 
FP= 1 • +F - 
FM* 1 . — F- 

AMC  = ( • 5/EXP 1  I  )*( AMCP*EXP1*FP+BMC*FM/EXPI  )- 

BMC  =  .5*EXPI  I*(AMCP*EXP1*FM4-BMC*FP/EXPI  )- 

REC-BEC/AEC- 

RMC*BMC/AMC- 

RMAG=. ABS. (REC+RMC  >- 

RRR ( KK )  =  • 25*RMAG*RM AG— 

thet=thet+dthet- 

WRITE  OUTPUT. 2* ( 01 ( LL ) . 02 (LL ) • 03 ( LL ) . 04 ( LL ) ) - 
WRITE  OUTPUT .2. ( ( RRR  t  1  )  . 1  =  1 • 1 • I .LE.KKK  )  >- 
R (LL ) =RRR ( 1 )- 

00  THR0UGH(S102).I=2.1.I*LE.KKK- 
PROV IDED (RRR (  I  ).L.R(LL) >.R(LL)*RRR< I  )~ 

CONT INUE- 
R 1 =R (  1  )  — 

R2=R(2)- 
R3=R (31  — 

R4*R ( 4 )- 
R5*R ( 5  )  — 

0  = . ABS*  (R2-R1  l  +  .ABS.  (R3-R1  M-.ABS.  (R4-RI  l+.ABS.  CR5-R1  ) 

RR(MM,-R1- 

WRITE  OUTPUT. 2- 

WRITE  OUTPUT. 2.  (01  (1  )  ♦  02  ( 1  >.D3(I  >*0.4(1  ))- 
WRITE  OUTPUT. 2. (R1 . R2 .R3.R4.R5 )- 
WRITE  OUTPUT. 2- 

PROV IDEDlRl .L.RR<MM-l  >).L=L+I- 
PROV IDED (L.G.5) .TRANSFER  TO (SI  60  >- 
01 ( 1 >=D1 ( 1 )+DDD*(R2-Rl >/D- 
02 ( 1  )  =02 < 1 )+DDD*(R3-Rl  ) ZD- 
03  (  1  )  =03 ( 1 ) +000* ( R4-R 1  )/D- 
04(1  )  =04 ( I  >+DDD»(R5-Rl  >/D- 
WRITE  OUTPUT. 2- 
CONT INUE- 

END  PROGRAM (FIRST )- 
DATA 

9.6  1.0 

•001  .000 

.089520  .046316 

5.0 


9.6 

.001 

•097800  .0002 


Fig.  4.  Digital  computer  program  for  optimizing  the  thicknesses 
of  a  reactive  wall  with  four  lossy  layers,  (cont). 
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